Gravitational wave background in perfect fluid quantum cosmologies 
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We discuss the gravitational wave background produced by bouncing models based on a full 
quantum evolution of a universe filled with a perfect fluid. Using an ontological interpretation for 
the background wave function allows us to solve the mode equations for the tensorial perturbations, 
and we hnd the spectral index as a function of the fluid equation of state. 



PACS numbers: 



I. INTRODUCTION 



The theory of cosmological perturbations relies es- 
sentially on two assumptions, namely that the back- 
ground is described by pure classical General Relativity 
(GR), while the perturbations thereof stem from quan- 
tum fluctuations, although they are subsequently evolved 
classically. Quite apart from the computational useful- 
ness of this scheme, this state of affairs is rather in- 
complete, and one would expect instead a fully quan- 
tum treatment of both the background and the pertur- 
bations to be achievable. In fact, the overwhelming ma- 
jority of classical backgrounds possess an initial singu- 
larity at which the classical theory is expected to break 
down. In recent years, many quantum background cos- 
mological models have been proposed, which share the 
attractiveproperty of exhibiting neither singularities nor 
horizons 0, 0, , leading the universe evolution through 
a bouncing phase due to quantum effects, and a con- 
tracting phase from infinity before the bounce. These 
new features of the background introduce a new picture 
for the evolution of cosmological perturbations: vacuum 
initial conditions may now be imposed when the Universe 
was very big and almost flat, and effects due to the con- 
tracting and bouncing phases, which are not present in 
the standard background model, may change the subse- 
quent evolution of perturbations in the expanding phase. 
Hence, it is quite important to study the evolution of per- 
turbations in these quantum backgrounds. The aim of 
the present paper is to provide a step in this direction by 
considering tensor perturbations in quantum minisuper- 
space background solutions. Interpreting the quantum 
theory in an ontological way 0, Q allows to define quan- 
tum scale factor trajectories, which can then be used in 
the second order tensorial modes perturbation equations 
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as show in Ref. . 

Note that such models may be viewed as alternatives to 
the standard inflationary paradigm Most known al- 
ternatives l9lllCI t o inflation present a primordial bounc- 
ing phase [TJ • Note that such a phase can also be 
seen as a complementary ingredient necessary for a com- 
plete cosmological scenario to make actual sense, i.e. not 
to be plagued with a singularit y or to avoid fac- 
ing any trans-Planckian problem [14| if, for instance, the 
bounce occurs at a scale such that all relevant cosmolog- 
ical scales now never went trans-Planckian. The bounce 
phase has recently been the subject of a lot of attention, 
in particular in view of the fact that in many instances, 
it was found to have the ability to modify the primor- 
dial spectrum of scalar perturbations, thus paving the 
way to confront them to the observational data In 
the case of bounces in quantum cosmological models, al- 
though the evolution equations for the perturbations may 
be constructed |16(, they are rather complicated due to 
the fact that the background does not satisfy classical 
Einstein's equations. Hence, all works in this area had 
to rely on a semiclassical approximation. 

In this paper, we calculate the gravitational wave back- 
ground spectrum produced at the bounce transition when 
this phase is described by a perfect fluid and the the- 
ory is fully quantized: this is the first time such a cal- 
culation, not involving any semiclassical approximation 
is performed. The restriction to gravitational waves 
stems from the fact that the perturbation equations for 
this type of modes can be substantially simplified, even 
when the background is quantized Q. Scalar and vec- 
tor modes, however, exhibit technical difficulties which 
have not been solved yet, so that tensor modes are, for 
the time being, the only modes that can be studied in a 
completely quantum way. 

The paper is organized as follows. Recalling in Sec. ITT1 
how the ontological interpretation allows a simple sep- 
aration between the background and the perturbations, 
we explain in Sec. IIIII how the Bohmian trajectories for 
the scale factor are derived, and discuss their generality. 
Then, Sec. IIVI which is the core of this work, provides 
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the tensorial modes indices in the known cases. We end 
in Sec. V with conclusions and discussions. 



II. THE MODE EQUATION 

The action we shall begin with is that of GR with a 
perfect fluid, the latter being described by the formalism 
due to Schutz [l7). i.e. 

(1) 

where £ m = (8nG N /3) 1 / 2 is the Planck length in natural 
units (H = c = 1), p is the perfect fluid pressure whose 
density p is provided by the relation p = top, lo being 



a constant. The metric g in Eq. Q is of the Friedman- 
Lemaitre-Robertson- Walker (FLRW) type, whose line el- 
ement we choose to be given by 

d.s 2 = N 2 (r) dr 2 - a 2 hys (r) ( 7ij + toy ) faftia? , (2) 

i.e. we assume it is perturbed to first order and restrict 
attention to tensorial perturbations only, with w , . = 

and w l i = 0, indices being raised and lowered by means of 
the background metric 7y- of the spacelike hypersurfaces 
(the bar denotes a covariant derivative with respect to 
this metric); the lapse function N(t), once fixed, defines 
the gauge. 

After inserting Eq. J2J) into the action JIJ, and perform- 
ing Legendre and canonical transformations, the Hamil- 
tonian up to second order reads (see Ref. |jj for details) 



H = NH, 
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(3) 



which is nothing but the Hamiltonian of Ref. 0] expressed for a perfect fluid. In Eq. © and in what follows, we shall 
denote by K, the spatial curvature (K. = 0, ±1 for flat, open and closed space respectively) in order to avoid confusion 
with the wavenumber k below. The quantities P ,IP J , P T are the momenta canonically conjugate to the scale factor, 
the tensor perturbations, and to the fluid degree of freedom, respectively. These quantities have been redefined in 
order to be dimensionless. For instance, the physical scale factor a p h ys can be obtained from the adimensional a 
present in J2J through a p h ys = £ Pl a,/"\/V, where V is the comoving volume of the background spacelike hypersurfaces. 
This Hamiltonian, which is zero due to the constraint Hq ~ 0, yields the correct Einstein equations both at zeroth 
and first order in the perturbations, as can be checked explicitly. In order to obtain its expression, no assumption has 
been made about the background dynamics. 

In the quantum regime, this Hamiltonian can be substantially simplified through the implementation of the quantum 
canonical transformation generated by 



U = exp(iG q ) = exp 



12 



(3a( 



(4) 



where (3 a = \ (^P a a + aP J and Q = J d' 
classical variables, yielding, for a particular factor ordering of (see Ref. for details) 



x^^WijW 1 ^ are the self-adjoint operators associated with the corresponding 
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(5) 



As we are here also quantizing the background, the quantization procedure is now to impose i?o^(a, Wy) = 0. The 
Wheeler-DeWitt equation in this case reads 



.5* 
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3u + l 



d J x 
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ij\k 



7 1 / 2 SwijSw^ 



24 



12 



*, (6) 
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where we have chosen T as the time variable, which is equivalent to impose the time gauge N = a 3uJ . Note that such 
a choice is possible in the case at hand because we are considering a a perfect fluid, for which one can use the variable 
which describes the fluid as a clock [HI El • 

Now, if one uses an ontological interpretation of quantum mechanics like the one suggested by dc Broglic 
and Bohm p|, and makes the separation ansatz for the wave functional ^>[a,Wij ,T] = <p(a,T)ip[a,Wij,T], with 
ip[a, Wij,T] — ipi[wij, T] J daip~ 2 (a, T) + ip2 [wi j , T] , then Eq. © can be split into two, yielding 

.dip a 3 "- 1 d 2 (p 



and 



i W = / d x 



-6- 



dT 

,30,-1) S 2 



da 2 



3w+l 
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ij\k 
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24 
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(7) 



(8) 



Using the Bohm interpretation, Eq. @ can now be solved as in Ref. yielding a Bohmian quantum trajectory 

a(T), which in turn can be used in Eq. ijHJ). Indeed, since one can view a(T) as a function of T, it is possible to 
implement the canonical transformation generated by 



U = exp 



d 3 i7 1/2 — ^ — 
2a 



exp < i 



d 3 x 



In 
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(9) 



where, as a(T) is a given quantum trajectory coming from Eq. I|7|. Eq. J^Jl must be viewed as the generator of a time 
dependent canonical transformation to Eq. (jHJ). It yields 



l ^ = I d 3 ^ 
dT 



,3w- 1 



1/2 ^ij\k^ 



ij\k 



27V 2 6^ij6^ \ 
Through the redefinition of time a 3w_1 dT = d?7, we obtain 



di] 



d 3 x 



27V 2 
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ij\k 
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(10) 



(11) 



This is the most simple form of the Schrodinger equation which governs tensor perturbations for a quantum minisu- 
perspace model with fluid matter source. 

The equation for the modes = to^/a which can be derived from Eq. I|ll|) reads (for that point on, the fc-index 
will be omitted) 



//' + r + 2K. Wt = 0, 



(12) 



which has the same form as the one obtained for classical backgrounds (see Ref. 0), with the important difference 
that the function a(r)) is no longer the classical solution for the scale factor, but the quantum Bohmian solution. In 
this way, we can proceed with the usual analysis, with the quantum Bohmian solution a (77) coming from Eq. Q acting 
as the new pump field. 



III. THE BACKGROUND BOHMIAN TRAJECTORIES 



In order to obtain the background quantum solutions, we choose the following factor ordering for the kinetic term 
of the background Schrodinger equation (from now on T = t): 



dp 



ft (3.-l)/2p na (3,-l)/2p (i 

d 



,3w+l, 



,(3u,-l)/2_ 



da 



,(3w-l)/2 



d_ 

da 



<p — Ka" ' ~tp 
p> — K-a 



if, 



(13) 



which is the factor ordering yielding a covariant Schrodinger equation under field redefinitions [20j. 

The quantum Bohmian trajectories are obtained from the solutions of Eq. I|13|) . There are two distinct situations: 
whether the spacelike hypersurfaces are flat or not. 



A. Flat spatial sections 

With the factor ordering chosen in Eq. (|13|) with K, = 0, 
we can change variables to x = f(l~ w ) 



-l a 3(l-o))/2 



obtaining the simple equation 

.dp _ 1 d 2 ip 
l ~di ~ \~d~x 2 ' 



(14) 
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Note that this is just the time reversed Schrodinger equa- 
tion for a one dimensional free particle constrained to the 
positive axis. As a and x are positive, solutions which 
have unitary evolution must satisfy the condition 



dx 



dip* 
dx 



(15) 



(see Ref. for details). 
We will choose the initial normalized wave function 



3 y/ 4 / x 2 



(16) 



where to is an arbitrary constant which determines the 
width of the Gaussian and hence the probability ampli- 
tude of initial scale factors. The Gaussian ipo satisfies 
condition (|15|) . It is a commonly used initial condition 
when the time gauge is fixed and one gets a Schrodinger 
equation of the type of Eq. l(T4"|l 0, 0, E3 , and even when 
the ti me g auge is not fixed when constructing wave pack- 
ets liEit 



Using the propagator procedure of Refs. 0-OI- we ob- 
tain the wave solution for all times in terms of a: 



ip(a,t) = 



St 


1/4 r 


exp 





-At^ 1 -^ 



9(i3 + t§)(l- W )2 



exp 



4ia 3 ( 1 ~ 



9(t 2 + t 2 )(l-^) 2 



arctan^ 



(17) 



Due to the chosen factor ordering, the probability den- 
sity p(a, t) has a non trivial measure and it is given 
by p(a,t) = a^ 1-3 ")/ 2 \ip(a, t)\ . Its continuity equation 
coming from Eq. I|13fl reads 



dp 
dt 



d_ 

da 



= 0, 



2 da 

which implies in the Bohm interpretation that 

. _ a^" 1 ) dS 
a ~~ 2 da' 



(18) 



(19) 



in accordance with the classical relations a = {a, H} = 
_L a (3u-i) Pa and Pa _ QS/da. 

Inserting the phase of l|17f) into Eq. (|19fl , we obtain the 
Bohmian quantum trajectory for the scale factor: 



B. Curved spatial sections 

In this case, only for io — ^ (radiation) are there an- 
alytic solutions available. Here, t = rj, and there is no 
factor ordering ambiguity in the kinetic term. Equation 
(0) [or ljT3|) ] reduces to the time reversed Schrodinger 
equation for harmonic or anharmonic oscillators. Now 
the condition for unitary evolution reads 



da 



dip* 
da 



(22) 



a=0 



the probability density p(a, t) is the trivial one, namely 
p(a,t) — \ip(a,t)\ 2 , satisfying the continuity equation 



a(t) = a Q 



(20) 



Note that this solution has no singularities and tends to 
the classical solution when t — > ±oo. Remember that we 
are in the gauge N 
time through 



Ndt = adrj 



drj 



{ait)} 3 ^ 1 dt. 



(21) 



The solution 

can be obtained for other initial wave 
functions (see Ref. 3). 



dt da\2da P > ' 



yielding the guidance relation 



IBS 
2 da' 



(23) 



(24) 



Given the same initial wave function as before, we ob- 
tain (see Ref. Q), 



1/4 
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x exp — i 



ICrj 2 ) \/~K.a 2 cos (^fjC-qj sin (^/lC.T]j 



/Ccos 2 (y^v) +sin 2 (v^J?) 



— arctan 
2 



sin ( 



7T 
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(25) 



(we change t to r) when a; = |). 

The Bohmian quantum scale factor obtained through the substitution of the phase of Eq. I|25|l into Eq. (|24|l reads 



a{rj) = a a 



(1 - /Cry 2 ) sin 2 (VtCr)) 

1 H 

K v 2 



= a 



cos 2 ( VlCr] 



sin 2 (VlCri] 



ICr, 2 



(26) 



For K. = and radiation, we can obtain the wave so- 
lution and Bohmian trajectories either by taking the re- 
spective limits from Eqs. (|17J) and l(2()|l or Eqs (|25(l and 
(|26(1 . The resulting Bohmian scale factor is 



a = do 




Note that, for the curved space section solutions be 
realistic, there must have a long epoch after the bounce 
when the scale factor recover its classical evolution and 
the curvature is negligible, i.e. when the scale factor in 
Eq. 1)26(1 can be approximated in some large interval of 
rj to a (ry) oc 7/ in order for the model to be compatible 
with standard nucleosynthesis and cosmological observa- 
tions. This can be accomplished if 770 <SC 1. It means 
that the initial wave function 1(1 6JI must be a very cen- 
tered Gaussian around zero. The flatness problem is then 
translated to the quantum cosmological language to the 
question: why an initial Gaussian wave function of the 
Universe is so centered around a null value for the scale 
factor? 



IV. TENSORIAL MODES PROPAGATION 



namely Eq. I|12l) . in the Bohmian picture with the scale 
factor assuming the form l|2(J|) or (|26|l . it is our goal in 
this section to solve this equation in order to obtain the 
gravitational wave power spectrum as predicted by such 
models. The first two subsections deal with the flat spa- 
tial section case for — | < u> < 1 . The final one treat the 



(27) 

curved spatial cases for u = 



A. Power spectrum for flat spatial section 



Our first task consists in going from the conformal time 
r\ to the more convenient time variable t stemming from 
the change (|21|) . With a dot indicating a derivative with 
respect to t, the mode potential reads 



a 



2(l-3tj) 



a . I a 

- + (1 - 3u) - 
a \a 



(28) 



Having obtained in the previous sections the propa- 
gation equation for the full quantum tensorial modes, 



and Eq. (|12|) transforms into 



V = 0, 



(29) 



in which we have defined v = o'' 1 3w V- Specializing to the flat K, = case and setting 

t k - -1 

i P iVtov , x=— and k~ ™^ - u " 3uJ -^ 

to 

r 



with ko — (toa^ 



we obtain 



d v 



dx 2 



2(3^-1) 

„2 (1 , 2\ 3 (l-c) 



k 2 (i + x 2 y 



1 

(1 + X 2 f 



which is in a useful form for the practical purpose of 
numerical resolution. We shall assume the usual vacuum 



v = 0. (30) 
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state initial condition for the modes, i.e. we set pj 

y/3L 



■ exp [-ik (77 - Thai)] 



(31) 



where ?yi n i is an arbitrary (and physically irrelevant as 
was checked numerically) constant conformal time, which 
we set to zero in what follows without loss of generality 
Figs. to 13 show the actual mode calculated numerically 
with Eq. (H1J). 

The power spectrum can now be defined as Q 



k 3 V h 



2k 3 \x 2 



leading to 



k 3 V h 



2k 3 - 2 I 



ao 



(32) 



(33) 



which, although in general being a time-dependent quan- 
tity, happens to be constant in the expanding phase for 
the time period we are interested in. Therefore, it suffices 
to solve Eq. (|30|l with the initial condition l|31[) to obtain 
the gravitational wave power spectrum we are seeking. 
This is how we obtained the figures. 



B. Piecewise approximation and matching in the 
flat spatial section case 

1. Asymptotic behaviors 

Given Eq. (|2U[1 and the relation l|21[l between the two 
time parameters, one has 



l-3w -T 

a 2-Tl 



3w - 1 3 
3(w-l)' 2'" 



(34) 



where 2F1 is a hypergeometric function and we have as- 
sumed a common origin for both times (i.e. r\ = for 
t = 0). This can be simplified by considering that one 
is mostly interested, either for setting initial conditions 
or for observing the resulting power spectrum, for times 
much larger than the typical bounce duration, i.e. for 
t ^> to. Then, one recovers the usual perfect fluid power 
law solution for the scale factor, allowing us to write 



ki] = ± 



3(1 



l + 3w 



CO)- 1+3 M 



(35) 



where the sign is to be determined by that of x. 

Eq. (|12|l with K. = has a potential for gravitational 



waves that is written V s , 



a" /a, and which can be 



expressed in terms of the x variable as 



Vgrav = A (1 + X 2 ) 3 <- 1 ) [(3W " 1) X 2 + 3 (w - 1)] 



(36) 



where the constant A is given by 



A= — 
9 



,l-3u 



(w - l)i 



For large values of ij, hence of x, and provided u> 7^ ^, 
one gets 



Vc 



grav 



2(3^+1) 

(3w - I) Ax s'— 1 ' . 



(37) 



which vanishes asymptotically for all cases of practical 
interest (—1/3 < u> < 1). 

The case uj = 1/3 is a very special and simple one: the 
time t is conformal time ?y, Eq. I|29l is identical to Eq. l|12fl . 
and a(rf) = ao[l + (ij/vo) 2 ] 1 ^ 2 - hi a different context, 
this same mode equation (with K. = 0) was treated in 
Ref. [2^, yielding an spectral index k 3 Vh cx k n T given 
by n T = 2. 



2. Matching points 

Let us determine the end of potential domination 
point, denoted by x M in what follows, i.e., the time at 
which k 2 = a"(x M )/a(x M ). This is 



9(1 

211 - 3u\ 







-1) 



2(l + 3u) 



(38) 



This point will also be used to match different solutions 
and thus propagate the mode through the bounce. 

It is interesting to note that the corresponding point 
for the evolution of v, namely the point obtained by an- 
nihilating the bracket in Eq. H30J) . called x cx n, is 



so that the ratio is 

■^cxit 



:).\ 



- 9 (i— 
cx fc 4 !+ 3 ^ <C 1, 



(39) 



which, for long wavelengths (k <C 1) is much less than 
one (recall that — 1/3 < uj < 1). Therefore, the modes we 
consider in the numerical solution have no time to start 
oscillating before reaching x m . 



3. Solutions 
Putting Eq. JjJZJ into E q- 021 and 

using <■'>.") I) , we arrive 
at the conclusion that, sufficiently far from the bounce, 
the perturbation mode satisfies 



2 , 2(3w-l) 



k 2 + 



(1 + Zuj)W 



/i = 0, 



whose solution is 

M = Vv \c 1 {k)H^ ) {ki 1 )+C2{k)Hl 2 \k7 1 ) 



(40) 



(41) 
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FIG. 1: Time evolution of the mode function v for the equation of state u> = 0.1 and wavenumber k 2 — 5 x 10 -8 . The full 
line is \v(x)\/ Vl + x^ — fi/a and thus provides directly the power spectrum. The symmetric curves are background functions: 
dashed is the conformal time potential a" /a as given by Eq. 1281 . dotted is the term proportional to k 2 in Eq. 1291 and full is 
(1 + x 2 )~ 2 . The horizontal thin straight line gives the value of k 2 used to compute the figures. The left panel shows the full 
time evolution which was computed. For x < 0, there are oscillations only in the real and imaginary parts of the mode, so the 
amplitude shown is a non oscillating function of time. It however acquires an oscillating piece after the bounce has taken place. 
The right panel is merely a zoom for smaller time scales also showing ±xm (the dotted vertical line) and x c ^ t (the full vertical 
line, indistinguishable on that scale with the axis). One clearly sees that even though the mode indeed starts oscillating, it 
does so on a timescale such that it is approximately constant all the way to xm- 




FIG. 2: Same as Fig. Q with ui — 0.3. The oscillations are visible on the left panel and the detailed view of the right panel 
emphasizes that there is no discontinuity on the mode. It also shows both matching points. 



with 



3(1 -u) 
2(3w + 1)' 



ci and C2 being two constants depending on the wave- 
length, f/^ 1 ' 2 ) being Hankel functions. 

This solution applies asymptotically, where one can 
impose initial conditions on the mode, as well as in 
the matching region for which VJr rav ~ k 2 , provided 
fc 2 C 1. Demanding the Bunch-Davies vacuum normal- 



3vr - 4f („+i) 



ization (|31|l then implies 

c\ = and c 2 = £ p 



The solution can also be expanded in powers of k 2 
according to the formal solution 0] 



/' 



A 2 (k) 



1-k 2 

v 



di] 



'/ 



a 2 (ry) 
.1 .1 a 



a 2 (77) dry 



"dry 2 pdfj , 2 ,,, 



, (42) 
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where A\ and A 2 are two constants depending only on 
the wavenumber k through the initial conditions. Ne- 
glecting the Oik 2 ) terms, for the expanding phase, the 
A 2 term is known as the decaying mode, and the power 
spectrum Q32JI can then be approximated accurately by 
a constant; this constant power spectrum is the one we 
are looking for. Although this form is particularly valid 
as long as k 2 <C a" /a, i.e. when the mode is below its 
potential, Eq. (|4*2")l should formally apply for all times. 
In the matching region, the 0(k 2 ) terms may give con- 
tributions to the amplitude, but they do not alter the 
fc-dependence of the power spectrum. 

For the solution (|2U[I , the leading order of the solution 
reads 



a 



A\ + A^oa^ ^ ar ctana; 

A . , 3(w-l) / 7T 1 

Ai - A 2 t a y 'I - + - 



(43) 



where in the last steps we have taken the limit x — > — oo 
and identified the leading orders in x, with A\ = Aiao — 
^aQ U '~ 2 toA2 and A 2 — — a^~ 2 t A 2 . Propagating this 
solution on the other side of the bounce, in the expanding 
epoch, yields the required power spectrum, i.e. the limit 
for x — > +oo, namely 



A 7T 3(0,-1), . 



— (ii - Tfl; 
an V 



(44) 



const " fl 

where we have taken only the constant part of the modes. 

4- Matching and spectrum 

In order to get the spectrum of gravitational waves pro- 
duced by our bouncing model, it suffices to match ji and 
p! at r] M , corresponding to x M given by l|38|) through (1351) . 



At this point, the mode function l|41|) and its derivative 
read 

C 

v(Vm) = 7/| and m'07m) = ^Vk, ( 45 ) 
where the constants C and D are given by 

C = C2Vk^H<?\kr hs ), (46) 

and 

D = C2(H?\kr> M ) 



+ y/k^ \Hl 2 \(k VM ) - Hl% (kr, M ) } , (47) 



with 



y/2|l - 3^| 

l + 3w 



This is also expressed as 



= —p and //(r? M ) = f)Vl, (48) 
Vk 

with C = ^/T a£ u ~ 1)/2 C and D = a ( ^ u)/2 D/^. 

Matching fi [Eqs. (|43|l and (|45|l ] and its derivative with 
respect to conformal time, namely fi' = a 1_3w t n ~ 1 d/i/dx, 
one obtains, to leading order 



Ax 
A 2 



3uj - 1 
3a (uj- 1) 
2 

3/3(1 -w) 



C-al u ~H aD 



~ 3(1- 

fc^+D, (49) 



- 3(a) -1) , % 

fc^+D , (50) 



with 



a 



9(1 -a;) 2 
2 II - 3^1 



-, -1/(1+3oj) 
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FIG. 4: Power spectra for different values of the state parameter ui. The dashed lines represent the approximation 1)5 lft with 
the power index given by Eq. (1521 . while the full lines are the spectra obtained by numerically solving Eq. 13011 . Shown on the 
right panel are solutions for u G {0.2, 0.3, 0.4, 0.5, 0.6}. 



and 



= 



9(1 



211 - 3w| 



The coefficients A\ and Ai contain each a power-law 
behaviors in k. Because ui < 1, the power in A2 [Eq. I|50(l ] 
is negative definite and that in A\ [Eq. 1)49)) ] is positive 
definite. Therefore, A2 is the dominant mode and gives 
the spectral index, while A\ provides the sub-dominant 
mode that happens, incidentally, to correspond to an un- 
altered propagation of the initial conditions. One then 
gets the spectrum Ipffl , and finally the spectral index n T 
writing 



k 6 V h oc k % 



and we get 



12w 

l + 'iuj 



(51) 



(52) 



Note that the limit u> — > | of Eq. (|5*2"|) gives the correct 
index for radiation (see Ref. [l^), although the calcula- 
tion, in this case, should not be valid; this is due to the 
expected continuity of the spectral index with the equa- 
tion of state. The spectrum as calculated numerically is 
plotted on FigQ]for various values of u together with the 
approximation 1)52(1. 

It is interesting to note that this result was also pre- 
sented in the semiclassical approximation (classical back- 
ground and quantum perturbations) in Ref. |23j |. In 
Ref. [23, the asymptotic behaviors both in the past and 
future infinities are two, possibly different, power laws for 
the contraction and expansion phases, whereas the type 
of bounces we studied here is restricted to equal asymp- 
totic behaviors, i.e., for \t\ ^> to- Since the potentials in 
the equations for \x are smooth and large compared to 



k around the bounce, it looks like the full quantum ef- 
fects and details of the bounce do not change significantly 
the main spectral features of the gravitational wave pro- 
duced. It would be interesting to verify if this results 
still holds for other bounces, e.g. those having different 
asymptotic behaviors and/or more complicated shapes 
of the potential for [i. In this last situation, and if the 
results of Ref. 12] apply, one would expect the actual 
spectra to be different. 



Power spectrum for curved spatial sections and 
radiation 



In this subsection we consider the power spectrum 
of tensor perturbations for quantum cosmological back- 
grounds with curved spatial sections. As mentioned in 
Sec. Ill, only in the radiation case one can obtain ana- 
lytic solutions for the quantum background. Hence, we 
will restrict ourselves to this fluid from now on. 

Inserting Eq. (|2l)|> into Eq. i|T2")l . and noting that k 2 = 
to 2 — 3/C, where m is an integer greater or equal to 3 for 
K, = 1, and a real number greater than zero for K, = — 1, 
we obtain 



% 2 



n -4- / j~fi — 

\ [r;2/Ccos 2 (V^?7) + sin 2 (V^7?)] 2 

The effective potential 



V eS = 



[i^K, cos 2 (VKrj) + sm 2 (VlCT])Y 



li = 0. 

(53) 

(54) 



has one maximum given by 1/t]q and goes to zero when 
r\ — > 00 for K. = — 1. It oscillates between 1/i]q and 775% 
which are respectively a maximum and a minimum pro- 
vided rjo < 1, when JC = +1. Indeed, as we have seen in 
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FIG. 5: Amplitude amplification of the gravitational modes in 
the K = — 1 curved case filled with radiation. For m<l, the 
amplitude behaves, as expected at least for the case tjo < 1, 
as a power law since then fi oc m~ 3 ^ 2 sin mq. 



Sec. IIIII in order for the background models to be real- 
istic, one must have rjo "C 1. Hence, the maxima of the 
effective potential are very high in both cases and the 
minima are very small in the K, = 1 case. Large wave- 
lengths (small to) will cross the effective potential and 
the perturbations will be amplified at each bounce. This 
induces an instability of the model because this enhance- 
ment happens an infinite number of times, and therefore, 
however small the initial perturbation might have been, 
there is a time at which the linear theory is no longer 
valid and the cosmological setup breaks down. 

For K. = — 1, the situation is very similar to the fiat 
case. The conformal time (ry M of potential crossing is 
given as the solution of the equation 



sinh(r? M ) = ±\ 



- 1 



mr) J ^1 + nl ' 



(55) 



which has a real solution provided 777770 < 1- Thus, the 
mode crosses the potential only in this case. This con- 
straint is however satisfied for the situations we are inter- 
ested in, namely, 770 <SC 1 and m <SC 1. One then obtains 



arcsinhl \ — 



(56) 



We have two limiting cases, namely 770 3> m and rjo <C m 
yielding, respectively, 



|rU«ln[2(VTO) 1/2 ] 



and 



(Vo/m) 



1/2 



(57) 



(58) 



The effect of the potential for /1 is to increase its am- 
plitude by a factor shown on Fig. as well as to mix 
the exponential terms. This can be easily seen by the 



following approximation. For rjo -C 1, the maximum of 
the effective potential at 77 = is very large while for 
77 3> t]q it behaves like ryp/sinh (77), which goes to zero 
when 770 — ^ 0. Hence the effective potential (15411 can 
be well approximated to a Dirac delta in this limit. Its 
integration reads 



V e «(ri)dr) 



1 





1 + 


f 1 H — 77 ) sinh 2 77 









l) sinh (277) 



1) + {rg + 1) cosh (2r?)] 



2?7o V ^0 



(59) 



which implies 



V c s(r])dr) = 1 



9 
Vo 



tan" 1 — ~ -?L + 0(770), 
Vo 2 Vo 

(60) 

where in the last step we have assumed that 770 <C 1. 

Hence, one can approximate the effective potential by a 
Dirac distribution [l2j ir5(i])/2riQ. The solution for 77 7^ 
is then simply Eq. Ij31|l with k substituted by to for 77 < 0, 
and n = Ae lmri + 73e _lm '' for 77 > 0. Demanding that fi 
be continuous across the potential and imposing Eq. 153(1 
then leads to another matching condition, namely 



M '(o+)- M '(o-) = — Mo). 



2770' 



(61) 



One then finds that 

A = - 
B = 



ITT 



V3L 



477oTO 3 / 2 ' 



ArjQirfl/ 2 V 777 



(62) 



and finally, in the long wavelength approximation for 
which m < 1, that /i cx to -3 / 2 sin 77777. This is exactly 
what is obtained numerically, as shown on Fig. Note 
also that when the curves reach the 777770 > 1 region, the 
amplitude is the initial one: the mode has not crossed 
the potential as explained above. 

Note incidentally that solution (|62|l is exatly the same 
as the one obtained in Ref. [23] and in the present paper 
for radiation with JC = 0, where we used the match- 
ing method. This is because also in this case one can 
approximate the potential to a Dirac distribution as 

V ~ Vo/v 4 when 77 ^> 770, which goes to zero in the 
limit 770 — > 0. However, for the parabolic scale factor 
0(77) = a [l + (77/770) 2 ] also treated in Ref. j^, which 
solutions are quite different from 1)62(1 . the potential is 

V = l/(77 2 + ?7q), whose limit for 77 ^> 770 is I/77 2 , inde- 
pendent of 770. Hence in this case, the effective potential 
cannot be approximated by a Dirac distribution, and the 
final spectrum is very different. We thus confirm that the 
power spectrum of perturbations through a bounce may 
depend significantly of the details of the bounce itself. 
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With the coefficients (|62ll . one can calculate the spec- 
trum 



m 6 V h 



2m i 



(63) 



for the two possible matching points (|57I58|) . yielding 
m 3 Vh oc m 3 ln 2 (m) and rr?Vh = m 2 . Note that, as ex- 
pected, the case 770 *C m yields the same spectrum as the 
flat case: the two scale factors are quite similar in that 
limit. 



V. CONCLUSION 

We have obtained the power spectrum of tensor pertur- 
bations in bouncing quantum cosmological models with 
a perfect fluid satisfying p — cup for flat spatial sections 
and — ~ < u> < 1, and for curved spatial sections with 
u) = |. For flat spatial sections, the spectral index for 
large wavelengths is n T — 12cj/(1 + 3uj). The positive 
curved spatial section model is unstable, while the neg- 
ative curved spatial section model amplifies the modes, 
changing the amplitude to a power index of n T 3 or 
ri T = 2, depending on the parameters. All cases lead to 
oscillations in the primordial spectrum. 

The most interesting case is the one of radiation, which 
is the best perfect fluid model for the early Universe (all 
particles are ultra-relativistic). For almost flat spatial 
sections we have n T ss 2, which is different from the 
predictions of inflation. Hence, this model can be po- 



tentially tested against inflation in future observations, 
specially concerning the polarization of the cosmic mi- 
crowave background, Planck observations, and gravita- 
tional wave detectors if we calculate the amplitude of 
these perturbations. 

The next step would be to calculate the spectrum of 
scalar perturbations of these models. The dynamical 
equations for scalar perturbations are not, however, as 
simple as Eq. Ijl2|l . The steps we have taken in Sec. (2) 
in order to arrive at Eq. 1)12(1 in the case where the back- 
ground is also quantized are not so simple in the case of 
scalar perturbations, specially due to the matter terms. 
This is work in progress. Attainment of the power spec- 
trum of scalar perturbations is crucial not only to test the 
model against WMAP observations, but also to calibrate 
and obtain the precise spectrum of tensor perturbations 
for possible comparisons with LIGO and VIRGO future 
data. 
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